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Q\ • The c < 1 and c > 1 matrix models are analyzed within large N renor- 

malization group, taking into account touching (or branching) interactions. 
The c < 1 modified matrix model with string exponent 7 > is naturally as- 
sociated with an unstable fixed point, separating the Liouville phase (7 < 0) 
from the branched polymer phase (7 = 1/2). It is argued that at c = 1 this 
multicritical fixed point and the Liouville fixed point coalesce, and that both 
fixed points disappear for c > 1. In this picture, the critical behavior of c > 1 
matrix models is generically that of branched polymers, but only within a 
scaling region which is exponentially small when c — ► 1. Large crossover ef- 
fects occur for c — 1 small enough, with a c ~ 1 pseudo scaling which explains 
numerical results. 
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I. INTRODUCTION 



Matrix models have proved to be useful tools to understand several important issues in 
string theory and quantum gravity Q. It is well known that large classes of random matrix 
models are equivalent, in the so-called planar limit (number of components iV — > oo), to 
discretized models of Euclidean 2-dimensional gravity, or equivalently to discrete non-critical 
bosonic strings. Critical points exist, which allow to take a continuum limit. Some of these 
models are exactly solvable, and their continuum limit is 2-D gravity coupled to matter with 
central charge c < 1, or equivalently bosonic strings in d = c + 1 < 2 (in a linear dilaton 
background). Explicit results are in agreement with other approaches, such as quantization 
in the conformal gauge (Liouville theory) or in the light-cone gauge ||, or topological 
gravity |J. The so-called double scaling limit 0] allows to sum over topologies of the 2-D 
world sheet, and to deal with non-perturbative issues in string theory. 

One of the still poorly understood issue is the nature of the so-called c > 1 barrier. No 
explicit solution is known for matrix models which describe c > 1 unitary matter coupled to 
2-D gravity. Moreover, the predictions of continuum theories for c < 1 become meaningless 
in the domain c > 1. For instance, the KPZ formula ]5],|3|,|4| for the string exponent 7 and for 
the scaling dimensions of matter operators lead to unphysical complex values when c > 1. 
This c > 1 barrier is not simply a technical problem, since it is related to the occurrence 
of tachyons for bosonic strings in d > 2. Alternative quantizations methods for strings in 
d > 2 lead only to partial results |§ , and Liouville theory allows only to construct consistent 
theories for c = 7, 13, 19 ||, which are nevertheless difficult to interpret in terms of physical 
(non-topological) strings. 

The nature of the c > 1 phase has been investigated by numerical simulations. The 
simplest case is dynamical triangulations in rf-dimensional space (2-D gravity coupled to d 
massless free bosons, with c = d > 1) [|iOHl4}| . For large values of d (typically d > 10), these 
studies point toward a branched polymer phase, with a string susceptibility exponent 7 = 1/2 
(besides "pathological" phases with highly singular curvature defects). This confirms semi- 
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rigorous and more heuristic arguments |T6|-[18| which shows that if 7 > the model is 



"generically" in a branched polymer phase with susceptibility exponent 7 = 1/2. However 
for smaller values of c (typically 1 < c < 4 ~ 5) numerical simulations point towards an 
intermediate phase with < 7 < 1/2, without discontinuity at c = 1. 

Another class of models consists in multiple Ising spins on dynamical triangulations 
^2] , [l~3|^3| - p5T,p8i^ll . At the order-disorder transition point, one recovers a c = n/2 theory 



(n being the number of spins). Here also for n > 2 but not too large 7 is found to increase 
smoothly with c from to 1/2, without discontinuity at n = 2. For large n, the ordered 
and disordered spin phases (with the pure gravity behavior 7 = —1/2) are separated by 
an intermediate branched polymer phase with 7 = 1/2, where the spins are already disor- 
dered [p^ - p8|j41|| • The transition between the ordered pure gravity phase and the disordered 
branched polymer phase is a branching transition with 7 = 1/3. A similar phenomenon 
occurs for g-states Potts models when q is large p9|j30[| , and is expected for 0(n) models 



when n is large |31| (but on flat 2-D space these models are non-critical at the order- disorder 
transition for q > 4 or n > 2). 

Analysis of c > models by series expansions give qualitatively similar results 
|T^,[T2|^,|S2| . So if for large values of c the behavior of all these models is generically 
(i.e. without fine tuning) that of branched polymers, for smaller values of c > 1 the situa- 
tion is still very unclear. This lead several authors to conjecture that for c > 1 but small 
there exists an intermediate behavior, with < 7 < 1/2, intermediate between the KPZ 
regime and the branched polymer regime |25|| , while it is also argued that very large finite 



size effects may occur as c — > 1. 

Another approach is the large N renormalization group (RG) introduced in |33| . The 



idea is that a change iV — > iV + 5N in the dimension of the matrix (the string coupling 
constant) can be absorbed into a change in the couplings gi —>■ gi + 5gi of the matrix model. 
This defines a renormalization group flow, with fixed points which describe the continuum 
limit. Calculations at lowest order in perturbation theory give a picture of the RG flow 
in agreement with exact results for c < 1 models. However the agreement is only very 



qualitative. The positions of the critical points and the values for the critical exponents are 
quite different from the exact results. For some c < 1 models one can use the equations 
of motion to write new non-linear RG equations which lead to the exact critical points 



and exponents [34|| , but this rely heavily on the solvability of these models, and cannot be 
extended to the interesting non solvable cases. In addition these equations are non-linear 
and not really in the Wilsonian spirit. 

In this article I want to show that the c > 1 barrier can be simply understood in this 
RG framework, if one takes into account the so-called touching (or branching) interactions 
in the matrix models. The scenario that I propose relies on already known results and on a 
few natural assumptions. It leads to a simple explanation of the characteristics of the c > 1 
barrier, and of the numerical results for c > 1, and might be of some relevance for higher 
dimensional models. 

Touching (or branching) interactions generate random surfaces which touch at isolated 
points, i.e. produce a microscoping "wormhole" connecting two world-sheets. In matrix 
models they correspond to products of traces of powers of the matrices. For instance, the 
simplest matrix model with such interactions, first introduced in [53 , has for action 



where $ is an N x N hermitean matrix. The Tr($ 4 ) term in the action generates the usual 
random surfaces made of squares, while the (Tr($ 2 )) 2 term allows two surfaces to be glued 
along an edge. It can be solved exactly, and its phase diagram is depicted on Fig.l . Starting 
from the x = 0, g = 1/12 BIPZ critical point there is a critical line which corresponds to 
pure 2-D gravity (c = 0), characterized by the string exponent 7 = —1/2. This line ends 
at the end-point x = 19/241, g = 361/5808, characterized by the exponent 7 = 1/3. Then 
is becomes a branched polymer critical line, with 7 = 1/2, which passes through x = 1/2, 
g = 0. These features persist for < c < 1 matrix models with touching interactions 
P^|^9]|. When the touching coupling constant x is switched on, the 2-D gravity critical 
line, characterized by the exponent 7 < given by the KPZ formula, ends at a end-point, 



characterized by the new positive exponent [|36|,fl0|,[41 
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and becomes for higher x a branched polymer critical line, with 7 = 1/2. There is an 
interesting conjecture [42,43|] that the continuum limit of the modified matrix model at 
the end point is still a Liouville theory, but with (some) positive gravitational dressings of 
operators replaced by the negative gravitational dressings. 



II. A SIMPLE RG CALCULATION FOR THE (TR($ 2 )) 2 MODEL. 



It is possible to study touching interactions with the RG approach of |33|. This is a 
natural idea, since (Tr) 2 terms are naturally generated by the RG flow at second order in 



perturbation theory. It is quite simple to apply the method of sect. 4 of |33| to the model 
(H), and to obtain the 1-loop RG flow equations. 
One rewrites the action ([]]) as 

S N m = Ntr[j(l-b)-g T ) + — - + N 2 p (3) 

with b an auxiliary variable, and p an additional trivial coupling constant, which normalizes 
the vacuum energy. Starting from an (N + 1) x (N + 1) matrix, one integrates over the last 
line v* and row v of the matrix, after rewriting the (v *v ) 2 term in the action in terms of the 
auxiliary variable a 

- ^-{v*v? - Ng (a{v*v) + ^) (4) 

and one obtains an effective action S c f[ = Sn + AS" for the remaining N x N matrix $. For 
large N the variation of action is 

(S 2 $ 4 Q \ (b 2 a 2 \ 

AS[^] = Tr f ^-(1 +b)-g— + Log(l +b + ga- g&) \ + N f - + g— + 2p\ (5) 

Assuming that 6~g~a;^(T~iwe may expand the second order in g, and replace a and 
b by their saddle point values in S[Q] + A5[$] 



a c = -(1 + g + gi) - ( x + 2g - xg - g 2 )^Tr (^pj 
1 /$ 2 \ 

We thus generate Tr($ 2 ), Tr(<£> 4 ) and (Tr($ 2 )) terms in AS. Rescaling $ in order to keep 
the Tr(<l> 2 ) term fixed, we find that the variation of the action amounts to a renormalization 
of g, x and p, which defines the 1-loop /3-functions 

&9 = -jjPb(g, x) ; P g = g-6g 2 - 2gx 
Ax = -jjP x (g, x) ; (3 X = 2x- 3x 2 - 6gx 

Ap = -i(1 P (g,x, P ) ; P P = -2p- \ + I + I + f + 3 -f (7) 

(8) 

Of course the fact that the action ([I]) remains closed under RG transformation is true only 
at 1-loop. The renormalization of p is unimportant at that stage, and is usually neglected. 

The corresponding RG flow in the (x, g) plane is depicted on fig. 2. The arrows goes from 
large iV to small N. Besides the Gaussian attractive fixed point (g,x) = (0,0), we recover 



the pure gravity fixed point of [Q at (g,x) = (1/6,0), which has one unstable direction 
, as expected. In addition, we find a purely repulsive fixed point at (g,x) = (0,2/3) and 
another critical fixed point at (g,x) = (—1/6, 1). The interesting conclusions that we draw 
are: (%) we indeed obtain a RG flow which mixes Tr($ 4 ) and (Tr($ 2 )) 2 terms; (ii) there 
are two critical points, i.e. two possible critical behavior for the model; and (in) the two 
corresponding critical lines are separated by a tricritical point. However, when compared 
with the exact phase diagram of Fig. 1, there are qualitative differences. We expect that 
the (0,2/3) fixed point corresponds to the branched polymer critical behavior, since the 
g = line has an enlarged U(N 2 ) symmetry, and is thus stable under the RG flow. In the 
first order calculation, it has two unstable directions and is thus multicritical, while from 
the exact result we find that the true (0,1/2) branched polymer critical point is critical, 
and should have only one unstable direction. One should nevertheless remember that the 
position of the critical points obtained by this first order RG calculation differs from the 



exact positions by ~ 50%, and that the error on the critical exponents, and on the stability 
of the fixed points, is very large. 



III. RG FOR C < 1 MODELS 

What lessons can be drawn from this simple calculation for c 7^ matrix models? We 
know from exact results that for c < 1 the phase diagram in the (g, x) plane (where g is the 
cosmological constant coupling and x the branching coupling) is generically similar to that 
of Fig. 1. If we assume that the large- N RG picture stays valid, the corresponding RG flow 
is schematically depicted on Fig. 3 . There are two critical fixed points with one unstable 
direction, the first one (A) corresponds to non-critical string , the second one (B) is on the 
g = line and corresponds to branched polymers. In between on the critical surface there 
is a multicritical point (C), with two unstable directions, which corresponds to the modified 
matrix model. 

The properties of the RG flow near the two fixed points A and C can be easily deduced 
from the exact results. Let us first consider the 2-D gravity fixed point A, and linearize 
the RG flow around A, by mapping the couplings (g,x,- ■ ■) into the renormalized couplings 
(i.e. the scaling fields) (t ,t x , ■ ■ •). Here t is the usual cosmological constant, and t x the 
renormalized branching coupling, other t's correspond to possible other couplings. In terms 
of the new couplings the singular part of the vacuum energy f[t] of the matrix model 



must obey the RG equation 
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f = (9) 



with the linearized /3-functions 

f3i[t] = AiU + Oit 2 ) (10) 

which defines the scaling dimension Aj of the fields t,. Let us first consider the model with 
no branchings (t x = 0). Since in the planar limit /[to] scales like A r2 to~ 7 , using KPZ scaling 
H (we restrict ourselves to the case of unitary matter) Aq is given by 



We can now easily deduce the dimension A x of t x . Indeed, taking one derivative with respect 
to t x amounts to insert one "wormhole" , i.e. to take two planar surfaces with one puncture 

r 1 2 

and to glue them at the punctures (see Fig. 3). This means that -M-f scales as 4^ -M-f 



and therefore that 



A, = 2(A + 1) = 2^-^ (12) 

One checks that (as long as c < 1) A x > 0, which means that t x is an irrelevant coupling 
and that wormholes can be neglected in the continuum limit. 

Let us now consider the multicritical fixed point C, associated with the modified matrix 
model. This fixed point has two relevant renormalized coupling, that we denote to an d t x . 



t is the renormalized cosmological constant, already considered in [42, 43 1 , but we must also 



introduce the renormalized branching coupling i x . In [42,43] it is shown that in the planar 
limit, if /[to, • • •] is the (singular part of the) vacuum energy of the ordinary matrix model 
in the continuum limit, the vacuum energy / for the modified matrix model is simply given 
by a Legendre transform 

/[t , ■ ■ •] = extremum w.r.t. t of (/[t , • • •] + t t ) (13) 

This general argument can easily be adapted to take into account i x , and gives 

f[i ,i x , ■■■] = extremum w.r.t. t of (/[t , • ■ ■] + t t + (t ) 2 i x ^ (14) 

(|j~3|) implies that the dimension of to is 



A„ = -2-A = -l-^/lz£ (15) 
(which is consistent with (|^)), while ([14]) implies that the dimension of i x is 



A x = -2- 2A = -2^^^ (16) 

One checks that as long as c < 1, A < A x < 0. to and i x are both relevant couplings, and 
the most relevant one is the renormalized cosmological constant to- 
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IV. RG FOR C = 1 MODELS 



We can now consider what happens at c = 1. Then 

A = A = -1 , A x = A x = (17) 

Assuming that there are no other fixed points in the vicinity of A and C, the only explanation 
of this behavior is that the two fixed points A and C merge at c = 1 into a single fixed point 
C, with one relevant coupling to, and one marginal coupling t x . Assuming that the branched 
polymer fixed point B is unchanged at c = 1, the corresponding RG flow is depicted on Fig. 4 
. This explains in a simple way the special features of the c = 1 theories. Firstly, the non- 
analyticity of the critical exponents as c — > 1 is simply reproduced by a regular variation of 
the RG functions as c — > 1. Considering (1 clS db small parameter, the RG functions 

m = ^ - ti (is) 

simply reproduce the c — > 1 singularity of A , A^, A and A z . Moreover this explains the 
logarithmic deviations to scaling at c = 1. Without fine-tuning of the branching interaction, 
i.e. in our language for t x < 0, the singular part of the free energy scales as / oc to/l°g(^o), 
and at the multicritical point, i.e. for t x = 0, / cx tg l°g(to) [§1)0]- This is simply reproduced 



with the /3-functions ([Tjf), if we take into account the p renormalization. Indeed, the term 
fip§pf = ftp m ^ ne equation (|9|) gives nothing but the additional right-hand-side r[g] 
discussed in which has therefore a standard RG interpretation. The correct RG equation 



is 



N m + ^w a + < 3 'w: + i 3 '£\f = <> (19) 

and the scaling at c = 1 is reproduced if /3 and (3 X are given by ( p!8|) and if we take 

P p [t] = csttl + ■■■ (20) 
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V. RG FOR C > 1 MODELS 



This scenario for the RG flows near c = 1 leads to a simple picture for the c > 1 phase. 
The pair of fixed points (A, C) simply disappear, and the RG flow is schematically depicted 
on Fig. 5. The only critical fixed point is then the branched polymer fixed point B, and 
the whole critical line in the (g, x) plane is attracted toward B. This means that even if one 
starts without touching interactions (x = 0) they are generated by the RG and drive the 
model in the branched polymer phase. This picture is in agreement with the results of the 
numerical simulations for large c, but also explains the situation for c small. For c — 1 > 
but not too large, the pair of fixed points become a complex conjugate pair of fixed point, 
not too far in the complex plane from the real critical surface. It is expected that there will 
be a strong slowing-down of the RG flow in this region of the critical surface close to the 
complex fixed points (depicted as the shaded area A on Fig. 5). In this region, the RG flow 
for c > 1 should look very much like the RG flow close to C for c = 1. This implies that if 
one considers the matrix model without touching interactions (x = 0), as the cosmological 
constant coupling g approach its critical value g c , there will be a large crossover region, 
up to some g' c smaller than (but close to) g c , where the RG trajectory takes-off from the 
critical surface in A, with an effective RG eigenvalue A c g- ~ — 1. In this domain, asymptotic 
scaling is not reached, and one expects that the free-energy singularity is characterized by 
an effective string susceptibility exponent 7 e ff = —2(1 + 1/A e ff) — 0. Only when g' c < g < g c 
the true asymptotic scaling is obtained, characterized by the branched polymer exponent 



If one uses the toy RG flow given by (|18|), one can estimate the size of the true scaling 
domain. One can still simplify the RG flow equation ([TJD into 



7 = 1/2- 



to 



(21) 



i = -e-t 2 



(22) 



(with e 



(c — l)/24 and where t denotes the derivative w.r.t the RG "time", i.e. the 
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logarithm of the rescaling factor) This approximation is valid for e small and t close to the 
c = 1 fixed point t = t x = 0. From (p2|), for e > the "time" needed to pass the t x — 
barrier (i.e. to go from the initial condition t x = — oo to the branched polymer fixed point 
at t x = +oo) is ~ tt/v^ and from (^1|) during this crossing the most relevant coupling t 
inflates by a factor exp(7r/y / e). Therefore if one start from t x ~ — oo, in order to reach 
the vicinity of the BP fixed point at t x ~ +oo one must start from a to ~ exp(— ir/^fe). 
This crude argument should nevertheless give the size of the true scaling domain, which is 
exponentially small when c — > 1 



This implies that the real branched polymer scaling should be unobservable in practice if 
c — 1 is not large enough, and that one observes the cross-over effective scaling with 7 ~ 0. 
This scenario should be quite robust, and is in agreement with the numerical results. 

VI. APPLICATION TO MULTIPLE ISING SPINS COUPLED TO 2D GRAVITY 

In the above section I have considered only the cosmological constant coupling g and 
the branching coupling x, assuming that it is enough to fine-tune g (i.e to make t = 0) 
to make the model critical. Thus the above discussion and the phase diagram of Fig. 5 
apply typically to the non-critical string in d > 1 space, that is to the Gaussian model 
where d massless free bosons are coupled to 2-D gravity. The situation is slightly more 
complicated for multicritical models, where more than one coupling have to be fine-tuned 
to obtain a critical theory. A typical example is the multiple Ising spins model, which has 
been extensively studied. I shall show that the R.G. flow picture advocated in the previous 
section applies as well to this case, and leads to a simple and natural explanation of the 
numerical results. 

Let me first consider the n = 1 case. In this discretized model, one considers a ferromag- 
netic Ising model on a dynamical triangulation (in practice the spins live on the faces of the 




(23) 
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surface). In the model without branching the model has two couplings, the usual string cou- 
pling g, and the spin temperature T (the symmetry is not broken explicitly by an external 
magnetic field). For fixed T there is a critical coupling g c {T) where the geometry becomes 
critical (infinite area). At this g c {T), at low temperature (0 < T < T c ) the spins are in an 
ordered phase and the fluctuations of the spins are not critical, so the fluctuations of the 
geometry are that of pure gravity (c = 0, 7 = —1/2). At high temperature (T c < T < 00) 
the spins are in a disordered phase, and the fluctuations of the metric are still that of pure 
gravity. At the critical point T c the fluctuations of the spins are critical, and it is well known 
that in the continuum limit one gets gravity coupled to the c = 1/2 minimal model (free 
fermions), with 7 = —1/3 f44| . 

It is easy to add touching interactions in this model, so that one gets a model with g, 
T and the branching coupling x as coupling constants. The phase diagram for such models 
was already studied in |45j . Let me give a schematic description of the critical surface: this 



means that g is adjusted to its critical value g c (T, x) where the fluctuations of the geometry 
are critical, and the structure of the critical surface is studied as a function of T and x. 
For T = the spins are frozen and for T = +00 they are decoupled from the geometry, so 
in both cases one recovers a phase diagram similar to that of the one matrix model as x 
increases. For x small the system is in the pure gravity phase 7 = —1/2, until it reaches the 
branching transition 7 = 1/3, and for large xit is in the branched polymer phase 7 = 1/2. 
Now if one considers the Ising critical point (T c ), as x is increased it spans a critical line 
T Ising (x), characterized by the c = 1/2 gravity behavior 7 = —1/3, until the branching 
transition multicritical point is reached. At this point the string exponent is 7i S mg = 1/4 
, as given by (0). For higher x the system is in the branched polymer phase, irrespective 
to the value of T, so the pure gravity + ordered spin phase (which occurs for small x and 
small T) is separated from the branched polymer phase by a critical line which starts from 
the T = branching transition point and ends at the branching+Ising point. Along that 
line the transition is just the branching transition, thus it is characterized by the branching 
string exponent 7 = 1/3- The same is true for the pure gravity + disordered spin phase 
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(which occurs for small x and hight T), which is separated from the branched polymer phase 
by another branching transition line, which ends also at the branching+Ising point. The 
corresponding phase diagram (on the critical surface) is depicted on Fig. 7. 

The corresponding RG flow and fixed points on the critical surface are easy to guess, 
and are depicted on Fig. 8. A is now the Ising fixed point, and C the branching+Ising fixed 
point. Besides the renormalized cosmological constant to, with dimension Ao = —6/7, there 
is another relevant coupling at A, namely the renormalized temperature t 2 , which is coupled 
to the energy operator e. Its dimension is A2 = —2/7. The renormalized branching coupling 
t x has dimension A x . = 2/7. The unstable fixed point C is at the end of the Ising critical 
line (t 2 = 0). It has three unstable directions, corresponding to the three relevant couplings 
to, £2 and i x , with respective dimensions A = —8/7, A 2 = —2/7, A x = —2/7. The T = 
and T = 00 planes are stable under the RG flow, thus one recovers a pure gravity f.p. A' 
and a branching f.p. C on the zero temperature line, and analogous f.p.'s A" and C" on 
the infinite temperature line. The RG flow lines going from C to C and C" correspond 
to ordinary branching transitions, and separate the ordered (O) and disordered (D) phases 
from the branched polymer (BP) phase. 

I now consider the multiple Ising spins model, by coupling n copies of the Ising model 
to 2D gravity, with the same spin temperature T. As long as n < 2 one expects from the 
KPZ-like arguments that at the order- disorder transition the spins are critical, and that the 
critical point correspond to gravity coupled to c = n/2 matter. The RG flow on the critical 
surface is still that of Fig. 8, with the dimensions of the couplings to, t x , to and i x given 
by (pi), (|T2"D, ( |15D and (|16"D, and with A 2 (n) = A 2 (n) < 0. Assuming as above that the 



RG /?-functions are regular as n — > 2, the f.p.'s A and C must coalesce into a single fixed 
point with two relevant directions (corresponding to to and t 2 ), and one marginal direction 
(corresponding to t x ). On the other hand as n — > 2 the spin fluctuations are not critical near 
the zero and infinite temperature lines, so one does not expect any qualitative change in the 
RG flows away from the critical n-Ising line. The corresponding RG flow is schematically 
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depicted on Fig. 9. 

For n > 2 but small the RG flow picture that one obtains is very interesting. The pair 
of fixed points A and C disappears in the complex plane, but the RG flow along the n = 2 
critical line is almost unchanged, until one reaches the shaded region where the n = 2 double 
fixed point was located. Therefore there is now a thin "funnel" separating the ordered (O) 
and disordered (D) phases. In this funnel one flows into the branched polymer (BP) phase. 
The transition lines separating these three phases are the RG flow trajectories going to 
the branching f.p.'s C and C". This implies that the transitions are standard branching 
transitions with exponent 7 = 1/3. 

The discretized models of n-Ising spins on dynamical triangulations have no microscopic 
branching interactions, and correspond to the horizontal x = line. The phase diagram 
structure that one obtains for n > 2 is precisely the structure which is obtained for large n 
28), or in the simplified model of ||41|| . The low temperature ordered phase is separated 



from the high temperature disordered phase by an intermediate branched polymer phase. 
The transitions are just ordinary branching transitions with 7 = 1/3. Moreover at these 
transitions one expects that the spin fluctuations are never critical. The argument that I 
used to estimate the size of the scaling domain in c > 1 models can be easily adapted to 
estimate the width of the funnel, that is the size AT (in spin temperature) of the intermediate 
branched polymer phase. One finds that it is exponentially small when n — > 2, 

* r-rn / constant \ 

and in addition the size of the scaling domain in the most relevant coupling g, for which 
the branched polymer scaling and the branching transitions are observable, still scales as 
(p3|). This explains why this scaling has not been observed in numerical simulation and 
series analysis for moderate values of n, and why in practice one measures effective scaling 
exponents < 7 e fr < 1/2 intermediate between the c = 1 scaling 7 = and the branched 
polymer behavior 7 = 1/2. 

Finally a similar picture is valid for q-states Potts models and for the 0(n) interacting 
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loops models. For q < 4 or n < 2 one expects a RG flow similar to that of Fig. 8. However 
at q = 4 or n = 2 one knows from the flat space results that the ordinary critical fixed 
point will coalesce with the ordinary tricritical fixed point, since in flat space the transition 
is first order for q > 4 or n > 2. So for the models coupled to gravity, I expect that as 
q — > 4 or n — > 2 the critical f.p. A will coalesce with an ordinary tricritical f.p. A tric and 
with the critical branching f.p. C, and that they all disappear for q > 4 or n > 2. The 
logarithmic corrections to scaling at c = 1 will be more complicated, but the RG flow for 
q > 4 or n > 2 will be qualitatively similar to that of Fig. 10. This explains the qualitative 
similarity between the behavior of the multiple Ising spins models, Potts models and 0(n) 
models at c > 1. However it seems that for all these models the fluctuations of the matter 
fields are never critical, and just drive some geometrical transitions of pure gravity. 



VII. CONCLUSIONS 

In conclusion, I have proposed a scenario for the c > 1 barrier in non-critical strings, 
based on the inclusion of touching interactions and on a RG analysis. This scenario is simple, 
generic, and in agreement with exact and numerical results. However, let me stress that it 
is still conjectural. In particular, further numerical studies are needed in order to test the 
scaling laws that it predicts for c ~ 1. 

Let me end with a few more general remarks on the role of touching interactions in 
discrete gravity. I have argued that it is necessary to include such interactions in 2-D gravity, 



since they are generated by the large iV RG procedure of [33[ for matrix models. In fact 
there is another equally valid reason to take them into account. In the dynamical lattice 
approach of 2-D gravity, there has been several (mostly numerical) attempts to perform 



real space renormalization p!^j24]|25|f46| -p0| . This amounts to replace a block of neighboring 
cells of the random lattice by a single, larger cell. A problem with such procedures is 
that they inevitably lead to touching points, since a small "bottleneck" can be replaced by a 
single vertex, thus generating wormhole-like configurations. Some real-space renormalization 
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procedures consist in removing such touching points, when they connect a small surface (a 



it is perhaps better to keep these configurations, and to include branching interactions from 
the beginning in these models. 

Another interesting question is to understand if touching interactions have a simply 
stringy interpretation. For instance, I do not know if the scaling dimension A x of the 
"wormhole operator" in the modified matrix model can be reproduced by that of a local 
operator in Liouville theory. 

These remarks may also apply to higher dimensional models, which have been used to 
discretize 3-D and 4-D Euclidean gravity. The numerical simulations on dynamical 3-D and 
4-D triangulations show a phase transition between a negative curvature phase and a positive 
curvature phase which bears similarities with the branched polymer phase of 2-D gravity. It 
is quite possible that touching interactions are needed to understand this transition. 



I am grateful to E. Brezin, I. Kostov and J. Zinn- Justin for useful discussions, and to I. 
Kostov for his comments and his careful reading of the manuscript. I also thank J. Ambj0rn 
for his comments on the preprint, which lead to several notable improvements. 
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FIGURES 
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FIG. 1. Exact phase diagram for the 1-matrix model with (Tr($ 2 )) 2 interactions. 
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FIG. 2. RG flow at first order in x and g for the (Tr(<I> 2 )) 2 model. 
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FIG. 3. Inserting a wormhole 
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FIG. 6. Schematic RG flow for c> 1 
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FIG. 7. Schematic phase diagram for n = 1 Ising model 
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FIG. 8. Schematic RG flow for n < 2 Ising models 
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FIG. 9. Schematic RG flow for n = 2 Ising models 
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FIG. 10. Schematic RG flow for n > 2 Ising models 
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